Anomalous Hall conductivity from the dipole mode of spin-orbit-coupled cold-atom 

systems 
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Motivated by recent experiments [Lin et at, Nature 417, 83 (2011)] that engineered spin-orbit 
coupling in ultra-cold mixtures of bosonic atoms, we study the dipole oscillation of trapped spin- 
orbit-coupled non-condensed Bose and Fermi gases. We find that different directions of oscillation 
are coupled by the spin-orbit interactions. The phase difference between oscillatory motion in 
orthogonal directions and the trapping frequencies of the modes are shown to be related to the 
anomalous Hall conductivity. Our results can be used to experimentally determine the anomalous 
Hall conductivity for cold-atom systems. 
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Introduction. — Transport phenomena play a cru- 
cial role in understanding and characterizing condensed- 
matter systems. Two of these phenomena, the Hall effect 
and the anomalous Hall effect (AHE) were both discov- 
ered in the late 19*'' century. That the Hall effect is due to 
the Lorentz force has been understood since those days. 
The AHE, a transverse voltage or current present in fer- 
romagnets in the absence of a magnetic field, is related to 
spin-orbit (SO) coupling and has proven much more chal- 
lenging to understand (for a review see Ref. [1]). Since 
SO coupling is responsible for the AHE, anomalous-Hall- 
like effects should also be present for particles that do 
not carry charge, and, indeed, such effects are observed 
for magnons |2j, phonons |3| |4j, and photons (51. Al- 
though these effects were observed using heat currents, 
and they should thus be called anomalous Righi-Leduc 
effects, their physical mechanism is similar to that of the 
AHE. 

In this Letter we consider the AHE in homogeneous 
and harmonically-trapped cold-atom systems. (The AHE 
was considered in cold-atom systems in the presence of an 
optical lattice in two dimensions by Dudarev et al. [5].) 
As the atoms are neutral, the AHE here refers to a mass 
current perpendicular to an applied force in the absence 
of a Coriolis force. (For cold-atom systems rotation and 
the resulting Coriolis force play the role of a magnetic 
field and the Lorentz force.) Our investigation is mo- 
tivated by the recent experiment by Lin et al. [7] who 
engineered spin-orbit coupling in a Bose-Einstein con- 
densate [SI [5] with lasers. This experiment is one of the 
latest achievements in studying phenomena known from 
solid-state physics in a cold-atom setting. Other exam- 
ples are the Mott-insulator-to-superfluid phase transition 
[lOj . Bardeen-Coopcr-Schrieffer superfluidity [TT], the 
Berezinskii-Kosterlitz-Thouless phase transition [T^l US] , 
and Anderson localization [TJ. Important features of 
cold-atom systems are that new regimes of physics (as 
compared to solid-state systems) can be explored, and 
their great amount of tunability. Furthermore, cold- 



atom systems are in principle disorder free and have a 
well-known microscopic description making it worthwhile 
to undertake a detailed comparison between theory and 
experiment, whereas in solid-state materials typically a 
multitude of effects play a role which makes modeling 
harder. 

In the case of the AHE, for example, the difficulty 
in understanding the effect lies in part in the inter- 
play between so-called intrinsic and extrinsic contribu- 
tions. Intrinsic contributions come from spin-orbit cou- 
pling effects in the bandstructure, whereas extrinsic con- 
tributions arise from disorder. A recent theoretical ad- 
vancement in the understanding of the AHE is the semi- 
classical description in terms of equations of motion for 
Bloch wavepackets [121 IS]. In this description the in- 
trinsic contribution to the AHE stems from anomalous- 
velocity contributions to these semi-classical equations 
of motion jl7| . In modern language, this anomalous 
velocity results from the Berry-phase curvature of the 
Bloch bands that in turn is determined by the topology of 
the band structure. The relation between bandstructure 
topology and the Hall conductivity was first emphasized 
by Thouless et al. [IS] , and has regained interest with the 
very recent discovery of topological insulators [19j . 

In a typical cold-atom experiment steady-state cur- 
rents are not readily created and transport coefficients 
can be measured only indirectly. In this Letter we show 
that the anomalous Hall conductivity can be obtained 
from the properties of the dipole oscillation of a cloud of 
spin-orbit coupled cold atoms that is trapped in an ex- 
ternal harmonic trapping potential. The dipole mode is 
a collective oscillation of the center-of-mass of the cloud. 
According to Kohn's theorem [5D] , the frequencies of the 
dipole oscillation are equal to the trap frequencies. The 
SO coupling, however, breaks the harmonic nature of the 
system and as a result Kohn's theorem for the dipole 
modes does not hold. We find that spin-orbit coupling 
modifies the oscillation frequencies and that different di- 
rections of oscillation are coupled by the spin-orbit in- 



teractions. The phase difference between oscillatory mo- 
tion in different directions and the mode frequencies turn 
out to be related to the anomalous Hall conductivity. 
This result can be used to experimentally determine the 
anomalous Hall conductivity for cold-atom systems. Be- 
low we detail the semi-classical Boltzmann approach on 
which our findings are based, determine the anomalous 
Hall conductivity for homogeneous non-condensed Bose 
and Fermi gases, and show how this conductivity can be 
obtained from the dipole oscillation of trapped atomic 
gases. 

Semi- classical equations of motion — We consider 
spin-1/2 atoms with mass m trapped in an external po- 
tential V^^ (x) in the presence of a generic spin-orbit cou- 
pling. The Hamiltonian is 



H 



P 

2m 



+ V''^x)-M{p)-T , 



(1) 



with p and x the momentum and position operators of 
the particles and r the vector of Pauli matrices. The 
last term describes the SO coupling, that for spin one- 
half particles is without loss of generality given in terms 
of a momentum-dependent effective magnetic field M. 

At the semi-classical level we consider in first instance 
the dynamics of the expectation values of the position 
X = (i), momentum, p = (p), and spin s = /i(t)/2 
degrees of freedom. We obtain the Heisenberg equations 
of motion 



X = 



P = 



p 2dM 
m h dp 

Qyex 



dx 
M 



S = S X 



(2) 
(3) 
(4) 



We proceed by assuming that the spin degree of freedom 
is much faster than the motion of the particles. Thus 
we let the spin follow the effective magnetic field M adi- 
abatically, and only allow for a small misalignment be- 
tween the spin and the effective magnetic field that is 
first order in time-derivatives of the orbital dynamics. 
This approach is essentially exact in the linear-response 
regime. Hence, we solve the equation for the spin degree 
of freedom Eq. Q up to first order in time-derivatives 
by sex J2i'^+\M\i'^^ f^)'Pi' withm(p(i)) the unit 
vector in the direction of M. For spins opposite to the 
field the result is — s. Insertion of the result for s into 
Eq. ^ gives [mill] 



Xk = 



de 



p,k 



dp 



kp X B{p) , 



(5) 



where the band index k distinguishes between atoms 
with spin parallel (-I-) or antiparallel (— ) to the 
field M. Furthermore, the dispersion is given by 
£p,fc = p^ /2m — k\M{p)\, and the vector field 



Bcip) = ?iEa,b6{..,,z} ^"'"^ {dm/dpa X dm / dp,) ■ m/2 
determines the anomalous velocity contribution. We re- 
fer to B as the Berry magnetic field. 

Boltzm,ann equation and anomalous Hall conductivity. 
— We proceed by calculating the anomalous Hall conduc- 
tivity for a homogeneous gas from the Boltzmann equa- 
tion for the distribution function fk{x,p,t) for atoms in 
band fc, that is given by 



dt 



Xk. 



dfk 
dx 



p- ^^ = , 
dp 



(6) 



where we ignored collisions as the intrinsic anomalous 
Hall conductivity does not depend on relaxation [1]. In 
the above, dsk and p are given by Eq. (Is]) and Eq. (Is]) 
respectively. We consider a steady-state situation with a 
constant applied force F — —dV^^/dx acting equally on 
atoms in both bands, and define the conductivity tensor 
"" by J = cr • i^, where j is the particle current density 
which is given hy j = Y,k I ljj^fk{p)xk- The anoma- 
lous Hall conductivity ctah is the off-diagonal component 
of this conductivity tensor. The solution of the Boltz- 
mann equation leads to the anomalous Hall conductivity 



TAH= Y^ 
ke{+-} 



{2T:hf 



kN{ep,k)B,{p) 



(7) 



where N{t) = [e('^-'')/'=«'^ ± l]-i with fceT the thermal 
energy and /i the chemical potential, is the Fermi-Dirac 
(+) or Bose-Einstein (— ) distribution function that ap- 
plies for fcrmions or bosons, respectively. The above ex- 
pression for the anomalous Hall conductivity is the intrin- 
sic contribution due to SO coupling effects in the band 
structure. In cold-atom atom systems there is, unless en- 
gineered [T4j , no disorder and thus extrinsic contributions 
are absent. 

So far we have considered a generic SO coupling. In 
order to make a connection with experiments we will now 
consider a Rashba-Dresselhaus [5T] form of the SO cou- 
pling so that the effective magnetic field reads 



M{p) 



/3 



a 



h^'^' 2 



(8) 



where a and /3 are the coupling constants for Rashha and 
Dresselhaus SO-coupling respectively, and A is a spin- 
splitting energy. We then find for the Berry magnetic 
field 



B{p) 



4(a2 _ /32)a;i2 



{A^h? + 4(a2 + /32)p2 _ i6a/3p^pj^)3/2 ' 



(9) 



Note that in the experiments by Lin et.al. [7] an equal 
amount of Rashba and Dresselhaus coupling was realized, 
i.e., a = ±/3. It follows that in this specific case B{p) — 
and that ctj^h = [22|. It is however experimentally 
straightforward to consider a more general SO coupling 
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FIG. 1. (Color online) The anomalous Hall conductivity for 
bosons (solid) and fermions (dashed) as a function of nA''. 
The lines correspond from top to bottom to a//3 = (oo, 2, 1.5) 
where o? + 13^ = mKjY? . The spin splitting energy A = 
0.2A;bT. 




a m/ft A 

(Color online) Dipole mode frequencies as given by 
The solid line is calculated for /? = 0, the dashed 
lines correspond to /3 = O.Sjn/ft^A. The spin splitting is 
A = 0.2fcBr. The number of particles Na = 1.8 x 10^ and 
temperature T — 200 nK. The inset shows the sin ip as & 
function of a 



FIG, 
Eq. (1121 



[7]. The anomalous Hall conductivity vanishes in the 
absence of a spin splitting A, in agreement with the fact 
that the AHE occurs in ferromagnets. 

In Fig. [T]we show results for the anomalous Hall con- 
ductivity of bosons and fermions. We only show the re- 
sults for Of > /3 since (TAH(±a, ±/3) = crAH(a,/3) and 
o'ah(Q!,/3) = — CTAH(/?,a)- The results shown for bosons 
are above the critical temperature for Bose-Einstein con- 
densation (this temperature depends on a,/3 and A). 
The anomalous Hall conductivity is independent of tem- 
perature in the degenerate (nA^ 3> 1) limit for fermions, 
with n the density and A = {2TTfi? /mk^TY/'^ the de 
Broglie wavelength, as expected. 

Collective Modes. — We now study the dipole oscilla- 
tion of an atomic cloud of Na atoms in an anisotropic 
harmonic trapping potential of the form V^^{x) = 
^ [uj'^{x^ + j/^) -|- uj^z^) , with uJr and uj^ the trapping fre- 
quencies. The dipole (^ = 1) oscillations are pure transla- 
tions of the cloud with no changes in its internal structure 
that are described by the equations of motion for the cen- 
ter of mass position Xq = ^ ^^ / da; / dp/ {2Tr h^fkXk 

and velocity Vq = ^ X)fe I '^^ I dp/(27r?i)'^/fcd;fc. Hence, 
we make the following ansatz for the distibution func- 
tion, fk{x,p,t) = nk{x - Xo{t),p - mvo{t)), where 
nk{x,p) = N{ep^k ~ f^{x)) is the Bose-Einstein or Fermi- 
Dirac distribution function in the local-density approx- 
imation, with fj,{x) = 11 — y™(a;). From the Boltz- 
mann equation we obtain the equations of motion for 
the center-of-mass coordinates. For small oscillations, 
we linearize these equations of motion resulting in 

xq^H ■va~'^V''{xQ)y.B- (10) 

mv^ = -VV^'^{xo) (11) 

where H is proportional to the Hessian matrix of the 



dispersion and B is the Berry magnetic field averaged 
over the trap, which are given by 



B^ 


. 1 V 


Hij = 


= 1 \- 



d-'a; 



d-^x 



d^p 

J2nhf 

d^p 

{2nh) 



knk{x,p)B{p); 
^nk{x,p)H,j^kip), 



with Hijk = md^epk/dpidpj. For the mode in the z- 
direction we obtain the result w = w^ as predicted by 
Kohn's theorem and expected since the SO coupling only 
affects the dynamics in the x — y plane. The modes in 
this plane have frequencies given by 



U}± 



■ UJr 



A± 



\ xy 



^xx-^yy ) 



A\ 



(12) 



with A = {Hj:x+Hyy+B'^m'^uj'^)/2. When there is no SO 
coupling, a, /3 = we find w = cj^ as predicted by Kohn's 
theorem. For nonzero SO coupling the double degeneracy 
of this mode is lifted. The eigenmodes of oscillation are 
given by a^o (0 ~ {xism{uj±t + ip±),X2sm{oj±t)) with 
smip± = mB,uj±/{H^y + m'^B'^^ujl)'^/^. In Fig. [2] we 
show the mode frequency 0;+ and angle ip+ as a function 
of Of, in the special case a = ±/? we have B^ = and find 
ifj. = as shown in the inset. Another special case occurs 
when the SO coupling is of the pure Rashba or Dressel- 
haus form. Then H^y vanishes which results in 1^ = 7r/2. 
The phase difference between the two different directions 
of oscillations is determined by the Berry magnetic field. 
We can relate the average of the Berry magnetic field over 
the trap to the anomalous Hall conductivity for a homo- 
geneous gas with a density equal to the central density no 
of the trapped cloud B^ ~ o'ah('^o)/'t-Oj where crAH(fio) 



is given by Eq. (l7|. This conductivity can therefore be 
experimentaUy determined by measuring the frequencies 
iO± or the phase differences (p± of the modes. 

Discussion and conclusions — We have studied the 
dipole oscillation of a trapped gas of spin-orbit coupled 
cold-atoms, and found that these oscillations can be used 
as an experimental probe for the anomalous Hall ef- 
fect. In the experiments by Lin ei.a/.[7] the SO coupling 
strength is q;,/3 ~ mA/H^ and the Zeeman spin-splitting 
is A ~ 0.2kBT, using a temperature of T ~ 200 nK. 
Taking values for a and /3 of this order, we find that the 
angle ip^ ~ 0.03 and that {uj± — ujr)/oJr — 10% which 
appear to be observable. 

Up to this point we have not considered collisions be- 
tween the atoms leading to damping of collective oscil- 
lations. The harmonic nature of our system is explicitly 
broken by the SO coupling leading to relaxation of the 
center-of-mass motion of the cloud (such relaxation is ab- 
sent when the SO coupling is zero and Kohn's theorem 
prevails). This can be described phenomenologically by 
adding a term —Vq/t on the right-hand-side of Eq. ( 11 ) 



which would lead to damping of the dipole modes but 
does not affect the anomalous Hall conductivity. The 
frequencies of the damped system are n± = uj± + ij± 
with the damping rate 7±, up to first order in 1/r, given 
by 



1 



7± 



H^ 



2ujI + 2H,, 



OJt. 



-B? 



-1 



We note that the relaxation time r can in principle be 
calculated from the Boltzmann equation but considering 
this, given the above remarks regarding its importance, 
is beyond the scope of the present paper. 

In the adiabatic approximation that leads to the semi- 
classical equations of motion, spin directions transverse 
to the magnetic field M{p) are taken into account ap- 
proximately as they give rise to the anomalous veloc- 
ity terms. One could go beyond this adiabatic approx- 
imation and consider the (2 x 2)-distribution function 
fcra'ip) that allows for all possible spin directions. We 
have checked, by solving the Boltzmann equation for this 
distribution function in the coUisionless limit, that our re- 
sults for the anomalous Hall conductivity and the phases 
ip± are not altered. 

Possible extensions of this work are to consider the 
partially Bose-Einstein condensed phase for bosons, and 
the situation without Zeeman spin splitting A. In the 
latter case the AHE is absent, but there will be a spin 
Hall effect [23l |2l] that can be probed via the spin-dipole 
mode. (The spin Hall effect for cold atoms was proposed 
by Zhu et al. [3S] for a cloud falling due to gravita- 
tion.) We also intend to investigate the effects of spin- 
orbit coupling on other collective modes, in particular 
the quadrupole oscillation. 
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